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Abstract 

O^ . We review recent investigations of the critical behavior of ferromag- 

^ ' netic q-state Potts models on a class of hierarchical lattices, with exchange 

1^} I interactions according to some deterministic but aperiodic substitution 

rules. The problem is formulated in terms of exact recursion relations on 
a suitable parameter space. The analysis of the fixed points of these equa- 
tions leads to a criterion to gauge the relevance of the aperiodic geometric 
fluctuations. For irrelevant fluctuations, the critical behavior remains un- 



o 
o . 

f^ changed with respect to the underlying uniform models. In the presence 

of relevant fluctuations, a non-trivial symmetric fixed point, associated 
j^ ' with the critical behavior of the uniform model, becomes fully unstable, 

and there appears a two-cycle of the recursion relations. A scaling anal- 
ysis, supported by direct numerical thermodynamical calculations, shows 
'■^ ' the existence of a novel critical universality class associated with relevant 

C , geometric fluctuations. 

o 

O ' The discovery of quasicrystals, and the design and investigation of magnetic 

superlattices, provided strong motivation for the analysis of spin models with 
aperiodic exchange interactions. Along the lines of the Harris ^ criterion to 
/\ ' gauge the influence of quenched disorder on the critical behavior of simple fer- 

^ , romagnetic systems. Luck [^ has proposed that sufficiently strong geometric 

fluctuations, associated with deterministic but aperiodic interactions, may also 
change the ferromagnetic critical behavior of the underlying uniform models. 
Luck's heuristic criterion has indeed been checked and confirmed in a number 
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of cases (including the quantum Ising chain 0] , and the two-dimensional Ising 
model g). 

In a series of recent publications |^, ^, |^, |8[ ^ |I^ , we considered ferromag- 
netic Ising and Potts models, with aperiodic exchange interactions according 
to a variety of two-letter substitution rules, on hierarchical Migdal-Kadanoff 
lattices |l^ . Taking advantage of the lattice structure, we were able to derive 
an exact expression for Luck's criterion, and to analyze the novel critical be- 
havior in the presence of relevant geometrical fluctuations. We now present an 
overview of these calculations. 

As an example, consider the period-doubling two-letter substitution rule, 

A -^ AB, (1) 

B -> AA, (2) 

so that the successive application of this inflation rule on letter A produces the 
sequence 

A^ AB-> ABAA-> ABAAABAB ^ .... (3) 

At each stage of this construction, the numbers Na and Nb, of letters A and 
B, can be obtained from the recursion relations 

with the substitution matrix 

M=(; I). (5) 

The eigenvalues of this matrix, Ai = 2, which is the period of the transforma- 
tion, and A2 = —1, govern most of the geometrical properties of the sequence. 
The total number of letters, at a large order n of these constructions, behaves 
asymptotically as A". The geometric fluctuations are of the order IA2I". It is 
then convenient to define the wandering exponent of the sequence, 

that expresses the asymptotic dependence of the fluctuations with the total 
number of letters, AA^ -- A^'^. 

The period-doubling sequence can be used to construct aperiodic spin models 
on hierarchical diamond-type lattices with 6 = Ai = 2 bonds and m branches. 
To go beyond the simple Ising model, and introduce an extra parameter to work 
with, we consider a g-state ferromagnetic Potts model, given by the Hamiltonian 




Figure 1: Three successive generations of the simple diamond lattice (b = m = 
2). Letters A and B indicate the exchange parameters (which are chosen ac- 
cording to the period-doubling sequence). 



where Ci — l,2,...,q, at all sites of the hierarchical lattice, Jij > 0, and the 
sum refers to nearest-neighbor pairs of sites (for q = 2, we recover the Ising 
model). The couplings can take only two values, J a and Js, associated with 
a sequence of letters produced by the period-doubhng substitution. In Fig. 1, 
we indicate some stages of the construction of this model on a simple diamond 
lattice (with & = Ai = 2 bonds and m — 2 branches). Note that the choice of 
the same interactions along the branches of the diamonds is indeed mimicking 
an aperiodic layered structure. 

We may also consider more general substitution rules and basic "diamonds" 
with m branches and b bonds along each branch, which leads to a lattice with 
the intrinsic or fractal dimension 



L> = In (mb) I In (6) . 



(8) 



In each one of these structures, aperiodicity may be implemented by a substi- 
tution rule of the form 



{A,B) -^ (^"iB*--"!, A"=^S''-"=) 



(9) 



with < ni, 722 < 6. These sequences are characterized by a substitution matrix 
with eigenvalues Ai = 6 and A2 = n2—n\^ which leads to the wandering exponent 



InjAal 
InAi 



In \n2 — n\\ 



(10) 



Now we decimate the internal degrees of freedom of the diamonds to write 
exact recursion relations for the reduced couplings. For the g-state Potts model 
on a lattice with 6 = 2 bonds per branch, and with interactions according to 
the period-doubling rule, given by Eqs. (nl) and (0), it is easy to write |q] 



/ xaXb + g - 1 
\xA + a;_B +gr - 2 



and 






9-1 



\2xA^q-2 



(11) 



(12) 



Xb 





Figure 2: Schematic representation of the parameter space and the flow dia- 
grams for irrelevant (a) and relevant (b) geometric fluctuations. The non-trivial 
symmetric fixed point of saddle-point character for irrelevant fluctuations (a) 
becomes fully unstable in the relevant case (b). 



where xa,b = exp {ql3JA,B), with (3 = l/ksT. 

In parameter space, besides the trivial fixed points associated with zero and 
infinite temperatures, there is a non-trivial symmetric fixed point, that governs 
the critical behavior of the underlying uniform model (see the sketch drawn in 
Fig. 2). The linearization of the recursion relations about this symmetric fixed 
point leads to the matrix form 



Ax' 



A - 



C [q) M^ 



AxA 
AxB 



(13) 



where C [q] is a g-dependent structure factor, and M^ is the transpose of the 
substitution matrix. The eigenvalues of this linear form may be written as 



and 



A, 



Ai = C(g)Ai = C ((?)&, 



C(q)\2^C{q)X',^C{q)b- 



(14) 



(15) 



where we have used the definition of w, given by Eq. (0). For irrelevant geo- 
metric fluctuations, in which case the critical behavior remains unchanged with 
respect to the underlying uniform model, it is required that Ai > 1, and IA2I < 1 
(see Fig. 2a). Therefore, we can write 

kl^C{q)b^by\ (16) 

with the usual thermal critical exponent associated with the uniform system, 

VT = TT^—, (17) 



2 - a„ 



where D, given by Eq. (||), is the fractal dimension of the lattice, and a„ is 
the critical exponent associated with the divergence of the specific heat of the 
uniform model [^ ^. From Eqs. (^ and (p^), we have the structure factor 

C{q)^h^^-\ (18) 

Inserting this expression into Eq. dl3) for the second eigenvalue, we have 

A2 = 6^^-i+", (19) 

which leads to an analog of Luck's criterion for a hierarchical lattice. As we 
have already mentioned, fluctuations are irrelevant for IA2I < 1. However, if 
IA2I > 1, which corresponds to a fully unstable symmetric fixed point, geometric 
fiuctuations are relevant (see Fig. 2b). Therefore, the critical behavior departs 
from the uniform case if we have 

|A2|>l^c^>l-— ^, (20) 

2 - a„ 

which is an exact criterion of relevance for the Migdal-Kadanoff hierarchical 
lattices. For the Potts model on the simple diamond lattice, with 6 = 2 bonds 
and m — 2 branches, this criterion can be written as 

q>g^ = 4 + 2^2 = 6.828427..., (21) 

where q^ is the same critical value for the relevance of quenched disorder ac- 
cording to a calculation by Derrida and Gardner p4| in the weak disorder limit 
(also, note that an > for a uniform Potts model on a simple diamond lattice 
with q > Qc). 

The recursion relations are so simple that it turns out to be relatively easy 
to carry out an analysis of the critical behavior under relevant geometric fiuctu- 
ations. It is not difficult to show that, for q > qc, in which case the non-trivial 
symmetric fixed point is fully unstable, there is a two-cycle that seems to be the 
candidate to define a novel class of critical behavior H . In terms of the second 
iterates of the recursion relations, of which each point of this two-cycle is a fixed 
point, it displays a saddle-point character, with stable and unstable manifolds. 
Let us use standard scaling arguments to analyze the novel critical behavior. In 
the thermodynamic limit, the reduced free energy per bond can be written in 
the scaling form |jl3, M 

fix)^g{x) + ^f{x"), (22) 

where g (x) is a regular function, x" is a second iterate of the recursion relations, 
b is the rescaling factor (number of bonds of the diamond hierarchical lattice), D 
is the fractal dimension, and we use the factor b^^ because we need two iterates 
to go back to the neighborhood of the initial point in parameter space. This 
equation has the asymptotic solution p5| 



where x* is one of the points of the two-cycle, Acyc is the largest eigenvalue of 
the linearization of the second iterate of the recursion relations about any one 
of the points of the cycle, P {z) is an arbitrary function of period 1, and the 
critical exponent a, associated with the specific heat, is given by 

a = 2 - 2- — = 2 - 2- — . 

in Acyc In Acyc 

All the calculations that we have performed, for a g-state Potts model on a vari- 
ety of Migdal-Kadanoff hierarchical lattices, with different values of 6 and m, and 
suitable substitution sequences, point out that the values of a from this scaling 
analysis are unequivocally different from the values of a^ for the corresponding 
underlying uniform models ||]. For example, for the period-doubling rule and 
the simple diamond lattice (6 = m = 2), with q — 7, the two-cycle is located at 
{xa,xb)i = (5.285..., 7.642...), and (x^, 2:5)2 = (6.697..., 4.750...), with eigen- 
values of the second iterate Ai = 3.933... and A2 = 0.985..., from which we have 
a = —0.0022..., which should be compared with a„ = 0.010.... With q = 25, 
the weakening of the second-order transition is even clearer: the two-cycle is lo- 
cated at {xa,xb)i = (6.942..., 234.34...), and {xa,xb)2 = (39.023..., 3.831...), 
with eigenvalues of the second iterate Ai — 4.243... and A2 — 0.343..., from 
which we have a = 0.0817..., to be compared with «„ = 0.404.... 

To test the validity of the scaling arguments, and the role of the two-cycle as 
the responsible for the new critical behavior, we have performed direct numerical 
analyses of the singularity of the free energy [g|. In these problems, it is well 
known that the reduced free energy can be expressed as an infinite series [^3|, |l5| . 
For the Potts model on the 6 = 2 diamond-type lattice, with the period-doubling 
rule, we have the series 



f {xa,xb) 



1 

^ (2m)' 

n=0 ^ ' 



iln(4"'-f4")-|-<Z-2)+iln(2x^ 



(24) 



where x)^b indicate the nth iterate of the recursion relations. As usual, we as- 
sume uniform convergence, and differentiate term by term to obtain the specific 
heat per bond. The critical temperature can be determined with high precision 
by making use of the existence of the paramagnetic fixed point at T = 0, cor- 
responding to xa,b = 00, which causes the apparent divergence of the series if 
summed without using any regularization trick. Fixing the parameters q and m, 
and also the strengths of J a and J^, the critical temperature thus calculated 
in fact locates xa,b in the attraction basin of the two-cycle. According to the 
expectations, for irrelevant geometric fluctuations, or for the uniform model, 
this method yields a critical temperature that locates the system in the stable 
manifold of the uniform fixed point. 

The singularity in the specific heat can be analyzed by fitting the numerical 
thermodynamic data to a function of the form A + B |(T — Tc) jT^ " over a 
somewhat arbitrary scaling region. For the irrelevant and uniform cases, we 



have always obtained very good fittings, in excellent agreement with the values 
of au predicted by the scaling theory around the uniform fixed point. Of course, 
these fitted values did not present any detectable sensitivity on the particular 
values of Ja and Jb- The Ising model (g = 2) on the simple b — m — 2 diamond 
lattice, with exchange interactions according to the period-doubling sequence, 
had been previously and independently analyzed by Andrade |lq|, with the 
same conclusions. For larger values of g, in which cases the scaling analysis 
predicts positive values of a (although, of course, smaller than the values of 
ttii), the fittings presented excellent agreement with the scaling predictions. 
For weaker singularities (mainly negative values of a), the fitted values were 
always somewhat bigger than the scaling predictions, with better agreement for 
increasing values of q. For b = m = 2, and g = 7, we obtained a = —0.005(4), to 
be compared with the scaling result a = —0.0022...; for q = 25, we obtained a = 
0.08(2), to be compared with a = 0.0817.... From the numerical calculations, 
we have seen an oscillatory behavior of the specific heat as a function of T. 
The period of oscillation is roughly given by the argument in Eq. (E3), with 
better agreement for increasing values of q. These log-periodic oscillations, 
which are well known phenomena associated with hierarchical structures, and 
indicate the need of a more refined scaling analysis, have also been detected in 
the calculations of Andrade ||l^, |l^ . 

In conclusion, we have shown the existence of an exact criterion to gauge the 
relevance of geometric fluctuations in a class of spin systems on Migdal-Kadanoff 
hierarchical lattices. For irrelevant fluctuations, the critical behavior remains 
unchanged with respect to the underlying uniform models. In the presence of 
relevant fluctuations, the non-trivial symmetric fixed point, which is associated 
with the critical behavior of the uniform model, becomes fully unstable, and 
there appears a two-cycle of the recursion relations in parameter space. A 
standard scaling analysis, supported by direct thermodynamical calculations, 
shows the existence of a novel universality class associated with the relevant 
geometric fluctuations. 
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